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The effects of electron-electron correlations on the low-energy electronic structure and their rela¬ 
tionship with unconventional superconductivity are central aspects in the research on the iron-based 
pnictide superconductors. Here we use soft X-ray angle-resolved photoemission spectroscopy (SX- 
ARPES) to study how electronic correlations evolve in different chemically substituted iron pnictides. 
We find that correlations are intrinsically related to the effective filling of the correlated orbitals, 
rather than to the filling obtained by valence counting. Combined density functional theory (DFT) 
and dynamical mean-field theory (DMFT) calculations capture these effects, reproducing the exper¬ 
imentally observed trend in the correlation strength. The occupation-driven trend in the electronic 
correlation reported in our work supports the recently proposed connection between cuprate and 
pnictides phase diagrams. 


Many recent studies have focused on the importance 
of the occupancy of the bands close to the Fermi level 
(Ep) in promoting the superconducting state of the iron 
based superconducting pnictides mm- Most interest¬ 
ingly, when the filling is expressed with respect to the 
half-filled 3d 5 bands, a striking similarity of the iron 
pnictide phase diagram to that of high-temperature su¬ 
perconducting cuprates has been pointed out HE], sug- 
gesting links between strongly correlated Mott insulat¬ 
ing behavior and unconventional superconductivity. In 
the iron pnictides, electronic Coulomb correlations are 
largely driven by Hund’s exchange coupling Hi] leading 
to a strong dependence even of the normal state proper¬ 
ties on doping [7] and possible orbital-selectivity in the 
strong coupling regime |8J. Theoretical models including 
the five-fold orbital manifold and Hund’s rule coupling in 
addition to local Hubbard interactions HUo] reproduce 
the main consequences of the change of filling induced 
by doping, in particular a regime with strongly enhanced 
effective masses when approaching the d 5 configuration, 
even at local Hubbard interactions much smaller than 
the critical one. However, while this picture works well 
in the case of substitution of Ba with K or of Fe with 
Co DU. it fails in case of isovalent substitution, e.g., for 
BaFe 2 As 2 -ccPcc 5 where the correlations change with x but 
the nominal filling of the correlated bands does not [2]. 


From the theoretical side, it has also been pointed out 
that nominal doping does not necessarily lead to a con¬ 
sistent change in band fillings [12]. Quite generally, the 
link between the correlation strength in the Fe-pnictides 
and the valence count remains elusive. 

In this letter, we emphasize the conceptual difference 
between the nominal valence count corresponding to the 
band filling and the effective orbital occupancies of the 
correlated states not only by doping but also by vari¬ 
ations of the hybridisations. Based on a combined ex¬ 
perimental and theoretical spectroscopic study, we show 
that the latter is a more reliable tuning parameter for 
electronic correlations. We perform SX-ARPES for a 
series of stoichiometric 122 pnictides, namely, LaFe 2 P 2 , 
CaFe 2 P 2 and BaFe 2 As 2 - The use of SX-ARPES, with 
its increased probing depth compared to ARPES per¬ 
formed with ultraviolet light as the excitation source 
(UV-ARPES), is essential for the determination of the 
bulk electronic structure of these materials since it has 
been demonstrated that they show a surface state when 
measured in the UV range [HHE]. Moreover, we com¬ 
bine the SX-ARPES results with DFT+DMFT theoreti¬ 
cal calculations to show how the electronic correlation ef¬ 
fects are associated with the detailed electronic structure 
and filling in the studied series of stoichiometric 122 pnic¬ 
tides. In particular, when d-orbital occupations are con- 
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FIG. 1: DFT band structure and SX-ARPES intensity maps of LaFe 2 P 2 - (a). Band structure along high symmetry lines. 
The thickness of the coloured line is proportional to the corresponding character of the band, (b) Fermi surface plotted using 
XcrysDen visualization package [18] (c). FS map with his — 575 eV. The superimposed dashed lines are the DFT FS sheets, 
(d). FS map along cut 2 in (c), taken with his — 440 — 640 eV in steps of 2.5 eV. The superimposed dashed lines in (c) and 
(d) are the FS from DFT calculations for k\\ <0. (e), (f). ARPES spectra along cut 1 in (c) at his = 575 eV and 525 eV, 
respectively. The cut at his = 575 eV ( his = 525 eV) is close to the Z-Y (Y-M) direction, (g), (h). ARPES spectra along cut 2 
in (c) at k z = 21 and k z — 20, respectively. The DFT bands are renormalized by a factor (Wdft/W) arpes = 1.5. 


sidered rather than nominal chemical valences a clear link 
between the evolution of the strength of electronic cor¬ 
relations and filling emerges. Indeed, contrarily to naive 
valence counting arguments, we find that the substitution 
of La with Ca does not change the d-orbital occupation, 
whereas the seemingly isovalent substitution of As by P 
does. Furthermore, we identify the source of the latter 
change in the occupation, namely the change in bond¬ 
ing (BB) antibonding splitting of the pnictogen atoms 
due to the isovalent substitution. Our results reveal that 
the changes in the correlation strength are mainly driven 
by unexpected changes in the orbital-resolved occupancy 
of the bands close to E Fj suggesting that such changes 
may be fundamental in determining the superconducting 
properties of a given compound. 

The SX-ARPES experiments were performed at the 
Advanced Resonant Spectroscopies (ADRESS) beamline 
at the Swiss Light Source (SLS). The experimental ge¬ 
ometry is described in m • Data were collected using 
circularly-polarized light with an overall energy resolu¬ 
tion of 50-80 meV. The samples were cleaved in situ at 11 
K and measured in a vacuum always better than 5 x 10 “ 11 
mbar. The k values are determined by taking into ac¬ 
count the photon momentum and are expressed in units 
of ( 27 r/a, 2tt /a, 2tt/c). The k z values were extracted by 
using the free-electron final-state approximation with an 
inner potential Vo = 13 for LaFe 2 P 2 and CaFe 2 P 2 , and 
Vo = 14 for BaFe 2 As 2 - All the intensity maps shown in 
the main text are obtained by integrating the ARPES 
spectral weight in an energy window of E F A 10 meV. 

Figure [lja) shows the band structure of LaFe 2 P 2 ob¬ 
tained using the Wien2k software package [17! . The lat¬ 
tice constants used in the calculations are listed in Table 
[T] The calculated Fermi surface (FS) [Fig. [ljb)] is similar 
to that of LaRu 2 P 2 presented in m with the exception 
that there is an additional band (7 in Fig. [l]) that crosses 
E f along the Y-X direction in LaFe 2 P 2 and forms a small 
hole pocket at the zone corner. 


In Figs. [ljc) and [ljd) we plot the ARPES spectral 
weight mapping at E F near the (k x , £ 7 , 21) and in the 
(£ 7 , 2—£ 7 , k z ) plane. The k z values were extracted by us¬ 
ing the free-electron final-state approximation (FEFSA) 
[23] . The overlaid dashed lines show the FS from DFT 
calculations for k x < 0 (The FS for k x > 0 is obtained 
by reflection with respect to the k x = 0 axis). The k z 
variation as a function of (k x , k y ) in the ARPES mea¬ 
surements with fixed photon energy has been taken into 
account according to the FEFSA. The observed FS is in 
good agreement with DFT calculations and it is highly 
three-dimensional, similarly to the case of LaRu 2 P 2 m- 
The ARPES intensity as a function of energy and k 
along cut 1 in Fig.[l][c) at his = 575 eV (k z ~ 21 ) and 525 
eV (k z ~ 20) are shown in Figs. []Je)-[lJf), respectively. 
Since the value of k z changes substantially along cut 2 in 
Fig. 0c) (A k z ~ 0.5), to trace the band dispersion along 
the Y-X and Z-A symmetry axes we made interpolations 
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FIG. 2: Curvature, EDCs and MDCs analysis for the Y—M di¬ 
rection. (a) Curvature intensity plot obtained from Fig. 0f). 
(b) MDCs and EDCs peaks supperimposed on the renormal¬ 
ized DFT bands (c)-(d) EDCs and MDCs from Fig. [ljf). 
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FIG. 3: DFT band structure and SX-ARPES intensity maps of CaFe 2 P 2 - (a) Band structure along high symmetry lines, (b) 
Three-dimensional FS sheets of different bands, (c) FS map with kv = 525 eV. (d) Intensity maps in the (k\\ — 2, 1 — k\\, k z ) 
plane, taken with kv — 440 — 640 eV in step of 2.5 eV. (e), (f) Spectra along the (k\\ — 2, 1 — k\\) direction with k z = 18 (Z-A 
direction) and k z = 17 (T-X direction), respectively. The superimposed dashed lines are the DFT bands renormalized by a 
factor ~ 1.5. 


between spectra at different photon energies to obtain 
the spectra with constant k z shown in Fig. Bg)- 0 h). 
The DFT band structure (dashed lines) are superimposed 
on the data. Contrary to the case of LaRu 2 P 2 , where 
DFT reproduces the measured band dispersion quite 
well, the overall agreement between the ARPES spec¬ 
tra and the calculated electronic structure of LaFe 2 P 2 
becomes reasonable only once the bands close to Ep are 
renormalized by a factor (Wdft/W) arpes = 1.5 d= 0.1 
(red dashed lines). Assuming a purely local and orbital- 
independent many-body self-energy to be at the origin of 
this renormalisation, we estimate a quasiparticle residue 
of Z arpes = 0-67 zb 0.05 [24]. All the bands close to the 
Ep are renormalized by this same factor, as evident from 
the curvature [25] . EDCs and MDCs analysis of Fig. |2j 

The decrease of the bare bandwidth, due to the change 
in the principal quantum number of the electrons close 
to the Ep , as well as the increase of Hubbard interac¬ 
tion U and Hund’s coupling J (see Table 0 upon isova- 
lent substitution of Ru with Fe significantly contributes 
to the increasing correlations in the system. However, 
the observed band renormalization factor in LaFe 2 P 2 is 
smaller than in the antiferromagnetic (AFM) compound 
(e.g. BaFe 2 As 2 ) and in the superconducting 122 Fe- 
based pnictides (e.g. Bai_ x K ;E Fe 2 As 2 ), indicating that 
the above changes induced by the Ru -A Fe substitution, 
are not the only mechanisms responsible for the increase 
of the strength of the electronic correlations observed in 
these systems. 

To investigate the effect of hole doping in the system, 
we performed a similar study for CaFe 2 P 2 , which has 
one more hole per unit cell compared to LaFe 2 P 2 and 
LaRu 2 P 2 - Similar to these compounds, the 3D FS of 
CaFe 2 P 2 is in good agreement with the DFT calculations 
[see Figs. ia)Bd)] , while the band structure has to be 
renormalized by a factor of (Wdft/W) arpes = 1.5 ±0.1 
[FigsU e)-(f)]. This result is quite surprising, since the 
bandwidth renormalization is not influenced by the for¬ 


mal decrease in the number of the electrons due to the La 
to Ca substitution, contrary to the expected behavior of 
a system moving toward half filling [29] . However this is 
explained once the multiband nature of these systems is 
taken into account. Comparison between Figs, [lja) and 
[3ja) shows that the substitution of La with Ca does not 
result in a rigid shift of all the bands. Rather, only one 
band [7 in Figure [IJa)] is fully pushed above Ep , while 
the others (with Fe d character) crossing Ep are qualita¬ 
tively unchanged. This is mainly due to the fact that the 
7 band is very sensitive to the La to Ca substitution since 
it is essentially a La (Ca) band hybridizing (strongly in 
the case of La) with Fe d and P p states. The filling of the 
d bands is thus barely changed even if the total number 
of electron in the unit cell is decreased by 1. We note here 
that the comparison between our finding and the mass 
enhancement measured by various quantum oscillation 
experiments (m/moFT in Table |TT|) shows that while the 
electron-phonon interactions are responsible for the main 
contribution to the mass enhancement in LaRu 2 P 2 , most 
likely causing its low-T c superconducting state, they are 
negligible in LaFe 2 P 2 and CaFe 2 P 2 , which are indeed not 


TABLE I: The lattice constants and internal atomic positions 
for LaRu 2 P 2 19] . LaFe 2 P 2 20 , CaFe 2 P 2 21 , and BaFe 2 As 2 

m ■ 


Sample 

a(A) 

c(A) 

Zp n 

LaRu2P2 

4.031 

10.675 

0.3593 

LaFe2P2 

3.841 

10.982 

0.3554 

CaFe2P2 

3.855 

9.985 

0.3643 

BaFe2As2 

3.957 

12.969 

0.354 
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FIG. 4: DFT band structure and SX-ARPES intensity maps 
of BaFe 2 As 2 . (a) Band structure along high symmetry lines, 
(b) Three-dimensional FS sheets of different bands, (c) FS 
map of BaFe 2 As 2 with his = 840 eV. (d) Intensity maps in 
the (fc||, 0, k z ) plane, taken with his — 720 — 930 eV in step 
of 10 eV. (e) Spectrum along the (kn % 1) direction with at 
his = 840 eV (k z ~ 30). The superimposed dashed lines are 
the DFT calculations renormalized by a factor ~ 2.1. 


superconducting. 

Having discussed the role of the La to Ca substitu¬ 
tion, it is now important to clarify the origin of the 
change in the bandwidth renormalization from CaFe2P2 
to BaFe2As2- The nominal filling, the screened Coulomb 
interaction U and the Hund’s coupling J [see Table [TT] 
are very similar in the two systems, suggesting similar¬ 
ities in the electronic correlations. However our mea¬ 
surements for BaFe2As2 [see Fig. [ 4 ], for which a sim¬ 
ilar analysis as above leads to identifying a bandwidth 
renormalization ~ 2.1 =t 0.2 ( Zarpes = 0-48 ± 0 . 05 ) in 
agreement with the existent ARPES measurements in the 
UV range m , indicate that there is a sizable increase 
in the bandwidth renormalization. Again the change in 
the bandwidth renormalization can be understood once 
the effective filling of the d bands is considered, instead 
of the nominal value. A first indication that the actual 
filling of the d bands is decreased in BaFe2As2 comes 
from comparing the DFT calculations. Focusing on the 
d xz ( d yz ) bands, the comparison between Figsj 3 la) and 
[ 4 ][a) reveals that the bands at -250 meV in CaFe2P2 are 
very similar to the ones at Ep in BaFe2As2, i.e., these 
bands appear shifted by 250 meV with respect to each 
others. The change in the band structure in this case 
is mainly due to the P to As substitution (similarly to 
the case of CaFe2P2 and CaFe2As2 m), which results 
in a change in the distance between the pnictide atoms 
Pn (Pn = P, As) of two adjacent Fe-Pn layers. In case 
of smaller distance (in CaFe2?2) the stronger bonding 
(BB) antibonding (AB) splitting of the Pn bands results 
in the AB bands being above Ep [in particular the band 



FIG. 5: DFT and DFT+DMFT calculation results, (a) DFT 
calculation in CaFe 2?2 as a function of the P-P distance, (b) 
Experimental and theoretical values of Z. (c) Orbital de¬ 
pendent occupation of Wannier d orbitals within DFT and 
DFT+DMFT. (d) Calculated values of U c rpa , J c rpa and 
Wdft • 


at 0.5 eV at the X point in Fig. §a) m- In case of 
larger distance, as in BaFe2As2, the splitting is smaller 
and the same band now crosses Ep. To compensate the 
increase of the carriers in this band, the other bands shift 
of towards Ep. A clear demonstration of this mechanism 
is shown in Fig. [ 5 ](a) , where the DFT calculations for 
CaFe2P2 are performed as a function of the P-P distance 
for a fixed P-Fe distance. 

To summarize, we have shown that while the weak 
strength of correlations in LaRu2P2 is mainly due to 
its larger bare bandwidth Wdft [see Table |n| and Fig. 
| 5 jd)] the evolution of correlations strength in LaFe2P2, 
CaFe2P2 and BaFe2 As2 is driven by the effective filling of 
the d orbitals crossing the Ep. The strength of the cor¬ 
relations remains constant in LaFe2P2 and CaFe2P2 due 
to the similar filling, and increases in BaFe2As2 when 
the filling decreases [see Fig. [ 5 jb)]. This result differs 
from what would be naively expected from the nomi¬ 
nal Fe 3d band filling, i.e., the same strength of cor¬ 
relations in the 3 d 6,5 LaRu2P2 and LaFe2P2 and big¬ 
ger (Wdft/W )arpes in CaFe2P2 and BaFe2As2, which 
are both nominally in the 3d 6 configuration. We un¬ 
derline that the observed occupation-driven trend in the 
correlations confirms the link between cuprate and pnic- 
tides phase diagrams. The average orbital doping from 
our DFT calculations are 0.29 and 0.28 for LaFe2P2 and 
CaFe2P2, respectively. These results would put them 
in the overdoped non-superconducting side of the unified 
phase diagram proposed in [T], where the correlations are 
not strong enough to allow the development of unconven¬ 
tional superconductivity and to display differentiation of 
bandwidth renormalization i.e., selective Mottness, for 
electrons with different band characters. 

The qualitative description above can be con¬ 
firmed by theoretical calculations using the combined 
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TABLE II: Critical temperature, ground state, bandwidth renormalisation with respect to DFT extracted from ARPES, mass 
enhancements reported in the literature, Hubbard U and Hund’s coupling J calculated from cRPA, bare bandwidth obtained 
from DFT , vertical pnictogen-pnictogen distance and nominal occupation of 3d-shell for the samples LaRu 2 P 2 , LaFe 2 P 2 , 
CaFe 2 P 2 , and BaFe 2 As 2 . 


sample 

To (K) 

Ground 

W oft 

m 

Ucrpa (eV) 

J cRPA (eV) 

W DFT (eV) 

^Pn-Pn C^-) 

nom 

n d 

state 

W ARPES 

171 DFT 








LaRu2P2 

4 ng 

low T c SC 

1 m 

2 (22 

1.89 

0.502 

9.4 

3.00 [22 

6.5 

LaFe2P2 


non-SC 

1.5 

1.8 US] 

2.035 

0.695 

7.55 

3.20 Em 

6.5 

CaFe2?2 


non-SC 

1.5 

1.5 [27] 

2.632 

0.723 

7.50 

2.71 EH 

6 

BaFe2As2 

142 [22] 

AFM 

2.1 

2-3 [28] 

2.572 

0.78 

7.4 

3.78 E2 
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DFT+DMFT method which gives access to the spec¬ 
tral properties of these compounds. The DMFT cal¬ 
culations were performed in the implementation of [ 32 ] , 
using the the hybridization expansion continuous-time 
quantum Monte Carlo algorithm [ 33 ] as implemented in 
the TRIQS toolkit [ 34 ]. Localized Wannier orbitals were 
built by truncating the expansion of the initial atomic- 
like orbitals to an energy window W, chosen here as 
W = [— 7 . 5 , 3 ] eV (see [ 32 ] for details). Calculations have 
been carried out at an inverse temperature of (3 = 40 
eV -1 , in the paramagnetic phases of the compounds. 
The calculations of the Hubbard interactions and Hund’s 
coupling (or more generally the Slater integrals Fo, F2, 
F4) have been performed using the constrained random 
phase approximation (cRPA) method in the implemen¬ 
tation of [ 35 ] . 


The values obtained for the interactions are shown in 
Fig. §b). The DMFT spectral functions for the vari¬ 
ous samples are presented in Fig. [6] (a)-(d). The DFT 
band structures, overlaid on the respective spectral func¬ 
tions are renormalized by a factor (Wdft/W) dmft- The 
resulting quasiparticle residues (denoted Zdmft ) agree 
well with the Z m values (m = z 2 ,x 2 — y 2 ,xy,xz,yz) di¬ 
rectly extracted from the self-energy on the Matsubara 
grid as Z m = {l - Im ^^A ^o]}" 1 [see Fig. ^b)]. 
The calculated behavior agrees well with the experimen¬ 
tal one; the correlations are negligible in LaRu2P2, in¬ 
crease in LaFe2P2, remain constant in CaFe2P2, and 
reach a maximum in BaFe2As2- Fig. [6] (e)-(h) shows 
the imaginary parts of the local many-body self-energy 
on the imaginary frequency axis, a quantity which in the 
Fermi liquid regime exhibits linear behavior at low ener¬ 
gies (and a slope related to the quasiparticle residues as 
above). These curves show that the DMFT calculations 
display the trend of correlation strength for our series of 
compounds in two ways: not only do the self-energies 
increase in overall value, with the steeper slopes for 
BaFe2As2 corresponding to the stronger quasi-particle 
renormalisations (smaller Z m ), but also does the ener¬ 
getic extent of the linear regime decrease. In BaFe2As2, 


coherent quasi-particles exist only on very low energy 
scales, while the other compounds display larger coher¬ 
ence scales. In the spectral functions, the signature of 
this effect is the overall quite substantial broadening 
observed in BaFe2As2, as compared to the other com¬ 
pounds. The small underestimation of ( Wdft/W ) in 
the DMFT calculations as compared to experiment might 
be due to the local nature of the DMFT method [36 or 
the absence of dynamical screening in the calculations 
HI EH [ 38 ]. We also note that there is a sizable change 
in U c rpa and J c rpa only upon the Ru to Fe substitu¬ 
tion and that, even in this case, it does not have a large 
influence on the correlations in our calculations. Indeed 
the calculated for the LaRu2P2 are slightly changed 
by an increase in U c rpa and J c rpa (values obtained us¬ 
ing the same U c rpa and J c rpa nsed for BaFe2As2 are 
Z z 2 = 0 . 90 , Z x 2_ y 2 = 0 . 91 , Z xy = 0 . 92 , Z xz , yz = 0 . 90 ). 
The evolution of the filling of the Wannier d orbitals, 
both before and after the DFT-DMFT calculations, is 
shown in Fig. [ 5 ][c). The filling of the orbitals is constant 
in LaFe2P2 and CaFe2P2 and decreases in BaFe2As2, in 
agreement with the qualitative conclusions drawn from 
the band structure calculations. 

We finally remark that a previous work has proposed 
CaFe2?2 as structural analogue of the collapsed tetrag¬ 
onal (CT) non-magnetic phase of CaFe2As2 [ 27 ] . How¬ 
ever the measured correlation strength observed here for 
CaFe2?2 is smaller than the value reported for the CT 
phase of CaFe2As2 [ 39 H 4 T] , indicating that there might 
be a fundamental difference between the two systems. 

In conclusion, we showed that the evolution of the elec¬ 
tronic correlations in the presented series of pnictides dif¬ 
fers from the expected behavior inferred from the nom¬ 
inal filling of the correlated orbitals. The experimental 
trend can be qualitatively understood based on the effec¬ 
tive filling of the correlated orbitals and quantitatively 
reproduced by our DMFT calculations. We also demon¬ 
strated that the observed occupation driven trend in the 
correlations supports the recently proposed link between 
the unconventional superconductivity in the cuprates and 
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FIG. 6: DFT+DMFT spectral functions and self energies, (a)-(d) Spectral functions calculated with the DFT+DMFT method 
along high symmetry lines for LaRu 2 P 2 , LaFe 2 P 2 , CaFe 2?2 and BaFe 2 As 2 , respectively. The DFT calculations superimposed 
are renormalized by the factors (^dmft) -1 , with Zdmft displayed in Fig. [H] (b). (e)-(h) Imaginary parts of the local many-body 
self-energy on the imaginary frequency axis. 


Fe-pnictides. Finally, the high sensitivity of the correla¬ 
tions to small changes in the filling opens new ways of 
tuning the strength of electronic correlations in transi¬ 
tion metal pnictide systems by systematic chemical sub¬ 
stitutions, pressure or constrains. 
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